The Effect of Thermal Fluctuations on Schulman Area Elasticity 
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We study the elastic properties of a two-dimensional fluctuating surface whose area density is 
allowed to deviate from its optimal (Schulman) value. The behavior of such a surface is determined 
by an interplay between the area-dependent elastic energy, the curvature elasticity, and the entropy. 
0^ ' We identify three different elastic regimes depending on the ratio Ap/As between the projected 

, (frame) and the saturated areas. We show that thermal fluctuations modify the elastic energy of 

stretched surfaces {Ap/As > 1), and dominate the elastic energy of compressed surfaces {Ap/As < !)• 
When Ap ~ As the elastic energy is not much affected by the fluctuations; the frame area at which 
the surface tension vanishes becomes smaller than As and the area elasticity modulus increases. 
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. I 

I. INTRODUCTION 

• Surfactant molecules are formed by chemically bonding two parts with different affinities to solvents - a hydrophilic 
O head group (usually ionic or polar), and a hydrophobic tail (usually composed of non-polar hydrocarbon groups). 
, ^ , . When these molecules are placed in a mixture of water and oil they tend to rest at the water-oil interface and decrease 
the interfacial energy substantially I'] . An interface that is saturated by surfactants (as in the case of some surfactant 
T-H , mixtures) has a nearly-vanishing surface tension^. Other interfacial systems which are frequently assumed to have 
^ ' zero surface tension are amphiphilic fluid bilayers [3| ■ Models ignoring surface tension effects (and relating the physical 
[ properties to the curvature energy only) have been very successful in describing many features of bilayer systems, 
such as the shapes of vesicles |3| and the entropic interactions between membranes 0. 

The problem with the notion of vanishing surface tension is that it can be understood in two different ways, 
. depending on the quantity in question. To demonstrate the two possible meanings of the term "surface tension" , we 
^f) ' consider a flat fluid membrane in an aqueous solution. Let us first assume that the total number Ng of amphiphilic 
, molecules forming the membrane is fixed, i.e., we ignore the exchange of molecules between the surface (membrane) 
' and the bulk (solution) 6]. This approximation may be relevant for highly insoluble amphiphiles, such as lipids, 
^ [ whose critical micelle concentration (cmc) is usually very low 7]. (Typical values of the cmc for bilayer- forming 
G ■ lipids are in the range 10~^ — 10~^° M Q). Lipid bilayers are also characterized by large "residence times" of the 
molecules within the aggregates |^ . The free energy of such a membrane depends only on the surface density of 
the amphiphiles S = N/A, where A is the area of the membrane. It achieves a minimum at a certain area density 
So which is determined by the intermolecular forces [lo| . Small deviations from the equilibrium density can be 



Q , described by a quadratic approximation of the Helmholtz free energy 

where Ag = NgJ^o is called the saturated area of the membrane, and Ka is the stretching-compression modulus. 
Typical values of Ka for lipid bilayers are in the range Ka ^ 10^ ergs/cm^ |^ For such values of Ka, any 
appreciable deviation of A from As would involve an energy cost much larger than any other relevant energy scale 
such as the thermal energy and the energy associated with curvature deformations. Consequently, it is frequently 
assumed that the total area is effectively constrained to As and, therefore, the surface pressure 



It is the vanishing of 11 for A — As which is usually described as the vanishing of the surface tension. 
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The ideas presented above are essentially due to Schulman [3, 1131 ; and so the saturated area is sometimes called 
Schulman area. However, a careful examination of Schulman's argument reveals that 11 (which we referred to as the 
surface pressure) does not comply with our usual understanding of the concept of surface tension. The quantity called 
surface tension is supposed to describe the free energy required to remove molecules from the bulk and create a unit 
area of the surface jlSj- It, therefore, cannot be discussed within a model that excludes the exchange of molecules 
between the membrane and the solution (like the one described above). This brings us to the second model: consider 
a membrane whose area density is fixed to Sq, but whose total number of molecules Ng may vary. The grand potential 
ilc of such a membrane can be described by 



where the chemical potential is a Lagrange multiplier that fixes the average number of surface molecules. From the 
above equation we learn that (a) the quantity 7, which truly deserves the name "surface tension" , is proportional to 
the chemical potential of the molecules and (b) that there is no special reason for 7 to be equal to zero. 

The differences between the above two examples can be viewed from a slightly different perspective: In both cases 
we deal with saturated surfaces. However, in the first example we deal with the response of the membrane to the 
variation of the area density of the molecules. The Helmholtz free energy flc depicted by Eq.(^ is nothing but the 
elastic free energy of the membrane, which takes a quadratic form near the equilibrium (reference) state. In the 
second case, on the other hand, we look at the process of creating more saturated surface. The addition of molecules 
to the surface is governed by the interfacial free energy fla described by Eq. , which is linear in A. Although we 
have just argued that 7 and not — H should be regarded as the surface tension of the system we will adopt, in what 
follows, a unified approach: We will define the surface tension of a flat surface as cr = dF/dA, where F (which we 
henceforth refer to as the free energy) is the relevant thermodynamic potential, i.e., F = Qc in the first case and 
F — flc in the second. (Note that flc and flc are not related by a Legendre transformation but rather describe two 
very distinct physical problems.) The differentiation with respect to A is taken while all the other state variables are 
held constant We get a — —H when we apply our definition of the surface tension to Eq.(Q, and cr = 7 when we 
apply it to Eq.Q. For a fluctuating surface we will use a generalized version of the above definition for a [see Eq.(0J, 
later in this section]. 

To fully describe the behavior of a thermally fiuctuating surface we need to include both the elastic free energy 
^ and the interfacial free energy © in our Hamiltonian. The above two models represent the two extreme cases of 
isolated elastic surfaces (first model) and incompressible surfaces (second model). There also exists a third model in 
which the surface is assumed to be both isolated and incompressible il5i 1I61 QJ^ . This last model has been successfully 
used in theoretical studies of bilayer vesicles 0]. As has been mentioned earlier, lipid membranes are usually considered 
as nearly-incompressible due to the separation of energy scales between the stretching and bending energies. As a 
result, the behavior of bilayer systems had been hardly studied in the framework of the first model, and very little is 
known 0, 0, about the effect of thermal fluctuations on the elastic free energy ^ (which have been described 
above for the case of a flat surface) and the corresponding surface tension . This should be compared to our very 
good understanding of the influence of temperature on the surface tension 7 associated with incompressible surfaces 
13, [ij, (20l I21L [2^ I23I |24| . The stretching elasticity contribution to the surface tension is usually described by the 
linear relation a = Ka{A — As)/As^ and is simply added as a correction to 7 |Tlll20i|. 

The aim of the present paper is to fill the gap existing in the literature and to study the effect of thermal fluctuations 
on Schulman area elasticity. Understanding of this subject should improve our analysis of the results of stretching 
experiments of water-insoluble lipid bilayers |llL l25l | . It is also important for understanding the results of computer 
simulations where frequently the number of molecules forming the surface is fixed |2^. Most of the early work on 
this subject, which is of a semiquantitative nature, is summarized in the review by de Gennes and Taupin 0. More 
recently |24| . there has been an attempt to introduce the stretching elasticity into the problem in a non-ad hoc way, 
i.e., without simply adding it as a correction to the surface tension 7. In Ref. |24j the behavior of a membrane with a 
fixed microscopic area was investigated, and the surface tension was associated with the optically visible area. Here 
we consider membrane which is elastic at the microscopic level and which is characterized by a finite compressibility. 

In this work we carry out a statistical mechanical analysis of the behavior of a nearly fiat fluid surface that spans 
a planar frame of a total area Ap (the "projected area") and consists of a fixed number of molecules Ng- We will 
calculate the free energy of a thermally fluctuating surface and extract the associated frame tension which is deflned 



na = ^lNs = ^iEqA = 7^, 



(3) 



by 



dF 
dAp' 



(4) 



The latter quantity gives the force per unit length experienced by the frame holding the surface. The frame tension 
coincides with another quantity (which is denoted by r) that has the dimensions of a surface tension (i.e., energy per 
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unit area) - the "g^ coefficient" |2l|, 123, |2g . Tfie surface tension r is tlie one that is measured in flickering experiments 
|25j | , where the dependence of the mean-square ampHtudes of the thermal fluctuations on the wavcnumber q is fitted 
to the formula 

where ks is the Boltzmann constant and T is the temperature. In order to investigate the elastic behavior of the 
system we first develop (in section ^ a general formalism which can be applied to any fluctuating surface whose 
energy can be written as the sum of an area-dependent term and a curvature dependent term. The validity of 
our formalism is demonstrated by applying it (in section ITTH) to determine the surface tension of an incompressible 
surfaces whose area-dependent energy is depicted by Eq.Q. The well known result |2l| for the temperature dependent 
surface tension is recovered. We then turn (in section llV)l to the main subject of the paper - the elastic behavior of 
compressible fluctuating surfaces. Our discussion of the subject is divided into three subsections: In subsection IIVAI 
we focus on the behavior of compressed surfaces with the projected area Ap smaller than the saturated (Schulman) 
area Ag. We will try to see the extent to which thermal fluctuations (which increase the total area compared to Ap) 
modify the Schulman elastic energy. In subsection IIVBI we discuss the behavior of stretched surfaces with Ap > Ag. 
We will examine whether the thermal fluctuations reduce or increase the (zero temperature) elastic surface tension. 
The behavior of surfaces with Ap ~ Ag is studied in subsection IIV Cl However, we must stress here that for very 
small values of the surface tension (which are encountered in part of this regime, close to the area where cr = 0) the 
elastic behavior is dominated by the interfacial free energy ^ and the associated surface tension 7 rather than by 
the elastic free energy The case when Ap = Ag was investigated by Brochard et al. Their interpretation of 
the surface tension had been reexamined 0] because it is based on the calculation of the partition function at a single 
projected area. In section we summarize the results and discuss the other factors that might influence the elastic 
behavior of the surface. 



II. THE FUNCTION Z{AA) 

We consider a nearly flat surface that spans a planar frame of a total area Ap. To describe the microscopic 
configurations of the surface we use the Monge gauge z = h{x, y), where h is the height of the surface above the frame 
reference plane . We assume that the Hamiltonian that describes the elastic energy of each configuration can be 
decomposed into the sum of two terms: 

n = ni{A[h])+ [ dsl-nH^. (6) 

J A 2 

The first term on the right hand side (r.h.s) of the above equation describes the dependence of the elastic energy on 
the total area A of the surface and has yet to be specified. The second term is the bending energy associated with 
the curvature. It describes the energy difference between flat and curved surfaces with the same total area. In Eq.Q 
the bending energy is expressed using the quadratic approximation of the Canham-Helfrich Hamiltonian |30l.l3lll3^. 
where dS is a surface area element, H = ci + C2 is the sum of local principle curvatures (the total curvature), and 
K is the associated bending modulus. For simplicity we will restrict our discussion in this paper to surfaces with no 
preferred (spontaneous) curvature, and to fluctuations which do not change the topology of the surface. For a nearly 
flat surface, i.e., when the derivatives of the height function with respect to x and y are small - hx,hy <^ 1, we have 

A~ Ap+ I dxdy \ {S/hf = Ap + AA, (7) 

and 

/ dSH^~ [ dxdy {V^hf , (8) 

Js JAp 

where the integral runs over the frame reference surface, and AA > [in Eq.Q] denotes the excess area of the 
surface due to the fluctuations in the normal z direction. 

Using Eqs.®-®, we write the partition function of the system as 
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where /3 = l/ksT. The above partition function can be also written in the following form 



Z = d (AA) exp [-/3Hi {Ap + AA)] 
Jo 

X Jv[h]6(^AA-J dxdy i(V/i)2^ exp -P^^ J dxdy^{V^hy 

POO 

= d (AA) exp [-iSHi {Ap + AA)] x Z (AA) , 

^0 



(10) 



where S is the Dirac delta function. The function Z (AA) can be identified as the partition function of a surface 
whose total area is constrained to the value A — Ap + AA. In other words, it reflects the probability density of the 
surface to have a total area A when the different configurations are weighted by their bending elasticity only. This 
probability density is determined by two opposite trends: Entropy favors strongly fluctuating configurations while 
the bending energy makes the configurations with moderate slopes more preferable. The contribution of the area- 
dependent term in the Hamiltonian ^ to the statistics is taken into account separately by the additional Boltzmann 
factor, exp [— /37ii {Ap + AA)], in the integrand in Ea. (|10|l . In order to calculate the function Z (AA) we use the 
Fourier space representation of the delta function 



-1 pioo 



(11) 



which when substituted in Ea. (|10(l yields 



1 niOG 

Z{AA) = / duj, 

V [h] exp \ 



dxdy 



(12) 



So far we have treated our surface as if it were a smooth continuous medium. One should not forget, however, that 
the Hamiltonian ® is derived from a more microscopic description. This means that the conformations of the surface 
cannot be defined below some microscopic length scale: In the frame tangent plane we must use a coarse-graining 
length I comparable to the size of the constituent molecules, so that an area element ("patch") of linear size I contains 
at least one (preferably a few) molecules. Our choice of / (which is somewhat arbitrary) defines a coarse-graining 
length scale \{l) for the normal displacements of the membrane patches, where \{l) — h{2n/mkBT)^/^ is the thermal 
de Broglie wavelength of a membrane patch, and m in the definition of A is the mass of the patch which is proportional 
to the (mean) number of molecules forming it. The bending modulus n is also scale-dependent [sSlMj. 

Details on length scale smaller than / are eliminated from our statistical mechanical treatment of the system by 
considering only those conformations of the surface h{x,y) — h{r) with wave vectors in the range |g| < A = 2^/tt/I: 



h{r) 




«1<A 



dqlXhqC 



xq-r 



(13) 



In the above equation hq is the amplitudes of the Fourier component corresponding to q p35| . The value of A had 
been chosen to set the number of modes included in the spectrum to be equal to the number of microscopic degrees 
of freedom N , i.e.. 



^ Ap 
P (27r)2 



ttA^ 



(14) 



For simplicity we use a circular, rather than a square, region in g-space. The sets of height functions h{r) in real 
space which can be generated by the modes included in the circular and square Brillouin zones are different only in 
the details on the microscopic length scale I. Therefore, both Brillouin zones are equally adequate in describing the 
macroscopic behavior of the system. With the Fourier modes representation, and using the identity 



df'e''-g~p)-r ^ (27r)2^(g-p), 



(15) 
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the excess area [Eq.(|7|l] is given by 



(16) 



where Aq = 2t\ j yJ~Ap ~ 0, and the mode 9 = has been excluded as it corresponds to a constant shift in h which 
does not contribute to the excess area. The second integral in Ea. (|12|l now reads 



V [/iq] exp ■ 



{2n) 



dq 



[u:q^+[3Kq^] 



Tracing over the surface profile hq in Ea. (|17|) . is straightforward, giving 



exp ■ 



Ar, 



{ujq^ + Pnq'^) 
2^ 



(17) 



(18) 



This result should be substituted back in Ea. p2|l . The first integral in Ea. p2|) can be evaluated in the thermodynamic 
limit using the method of steepest-descent. We find that (up to a constant factor) 



where 



satisfies 



Z{AA) = exp{/3G(AA)} 
~ exp < Lj si A A) A A 



ojs = /3k 



Ar, 



dq 



(w, (AA) g2 + pnq^) ;2^2 



- A^ exp (8ttPkAA/Ap 
expiSnPnAA/Ap) - 1 



2n 



dG 
duj. 



-0, 



i.e., solves the equation 



AA- 



dq 



(2^)2 2(a>,+/3«;g2) 



= 0. 



(19) 



(20) 



(21) 



(22) 



If we now substitute our expression for Z (A A) H19|l in Eq. 110|l and again use the saddle-point approximation, we 
obtain 



Z ~ exp{/3 [G (AA*) - Hi (Ap + AA*)]} , 
where AA* is the solution of the equation 

dHi dG dG dojs 



dAA d{AA) dujsd{AA) 



= 0, 



(23) 



(24) 



and where the function G (AA, utg (AA)) is defined by Ea. (|19|l . The last term in the above equation vanishes by virtue 
of Eq.lETl, while from Eq.lO we have dG/d (A A) = Us/P- We thus find that A A* is the solution of 



Us (AA) = p- 



dHi 



d{AA) 

The free energy of the surface is given by 

F = -fcsTln (Z) = Hi (Ap + AA*) - G (AA*) . 



(25) 



(26) 



The first term in the above expression is the area-dependent energetic contribution to the free energy. The second 
term can be regarded as the entropic part of the free energy, where the different configurations of the surface are 
weighted by their total bending elasticity. 
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III. INCOMPRESSIBLE SURFACES 



In order to evaluate the free energy of the surface we now need to specify the area-dependent part Tii of the 
Hamihonian (jSJ. Let us first consider the extensively studied example of an incompressible surface (see discussion 
in sectional. In this case, the area density of the surface is fixed, and any change in the total area due to thermal 
fluctuations must be matched by a change in the number of molecules that fixes the area per molecule. The fluctuations 
in the area are governed by the Hamiltonian [see Eq. ©] 



Hi = 7A = 7 (Ap + AA) 



(27) 



where (as discussed in sectional 7 is directly proportional to the chemical potential for the addition of molecules to 
the surface. From Eas. (|25|l and (|27(l we find that 



Using Eas. (|19|l . (jSEJ, and H28|l we obtain the expression for the free energy (consistent with |2ll |) 



(28) 



"/Ap + kBT 



An 



ksT P dq 



In 



(/37g" + /3«g^) PX^ 
2n 

2^ 



(29) 



which is correct to the lowest order in an expansion in (/3k)^^ — ksT/n. To find the surface tension, cr, we need 
to take the full derivative of F with respect to Ap. It is therefore important to examine the implicit dependence on 
Ap of the expression inside the braces in Ea. H29|l . One part of this expression that obviously depends on Ap is the 
lower limit of the integral Ao = 27r/ Ap. This dependence is associated with the logarithmic finite size correction to 
the free energy, whose origin is the fact the variation of the linear size of the frame {y/ Ap) leads to changes in the 
wave-numbers of the long wavelength modes. The contribution of this effect to a is negligible in the thermodynamic 
limit. A more subtle issue is the possible dependence on Ap of the microscopic length / and of the quantities A and A 
which are directly related to I. Recall that the patches have been identified as small sections of the surface containing, 
on average, a given number of (at least one) molecules. Since the area density of the surface is fixed, / can be derived 
from the constraint 



N 



2^4) 

An 



Const. 



where {A) is the mean total area of the surface. Using Ea. (|29|l and the relation (A) — dF/dj, we get |2C 



{A) =Ap\l 



Stt/Jk " 



■In 



7 -I- kA^ 
7 + kA2 



(30) 



(31) 



For 7^0 and in the limit — > 00 ( Aq = 27r/y^A^ 0), we readily conclude that criterion H3()|l is obeyed by setting 
Z to a constant value which does not depend on Ap. We thus find (ignoring the above mentioned finite size correction) 
in agreement with '2H that 



dF 
dAn 



F 



7 + 



dq 



(2^)'Ao 2 



In 



(/37g2 + pnq^) PX^ 
2^r 



(32) 



If 7 = we find / depending on Ap, what brings in another finite size correction to the above result (|32|l . 



IV. COMPRESSIBLE SURFACES 



We now turn to study a surface with a fixed number of molecules, whose area elasticity can be approximated by 
the harmonic form [see Eq.QJ] 



n,^^JA-A.f^^JAA-AA,)\ 



(33) 
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where AAg = Ag — Ap. The system under consideration is characterized by three energy scales: fcsT = k, and 

Kjitts, where = Ag/N is the Schulman area of the surface patches (the microscopic degrees of freedom). In what 
foUows we will assume that 



yJPKAag > /3k > 1. (34) 

These relations are obeyed by typical values of the elastic moduli and the Schulman (saturated) area per lipid of 
phospholipid bilayers. We are interested in calculating the free energy of such a surface as a function of the ratio 
Ap/Ag between the projected and the saturated areas. Combining Eas. (|2Uf) . (|25l) . and (|33|) we arrive to the following 
equation 

exp (87rpKA^/Ap) — 1 Ag 

whose solution AA* should be substituted in expression (|26|l for the free energy of the surface. Three regimes can be 
distinguished: 



A. Floppy surfaces 



Consider the case when Ap < Ag. The elastic energy caused by the mismatch between the frame area and the 
saturated area can be relieved by thermal fluctuations which store the extra area needed to bring the total area A 
close to Ag. Let A* be the projected area for which the solution of Ea. (|35|l coincides with Schulman area 



so that 



One can easily verify that 



A* - Ag 



ujg = 0. 



Ap — Ag 



= Ag 



In - 



(36) 
(37) 

(38) 



For Ap ~ A* the corrections to Eqs.(|2SIl ^"^^ can be expressed as a power series in the variable eA = Ap — A* 
(le^l < [max(ln A^, 87r/3«;)]~-^A*). The following relations can be derived: 



A* = Ag 



{2Tr)'^KAg 
KAAf 



In 



and 



(27r)2/3K 
A* 









+ 1 - 8TTf3K 







1 - 8tt(3k + O 



Ap 



£A_ 
A% 



(39) 



(40) 



which when used together with Eas. (|19|l . (|26(l . and (|33(l . yield (after some straightforward, but lengthy calculation) 
the following expression for the free energy 



F = Fo+Fi 



knT 



In 



(27r)2 
1 - Stt/^k 



dq 



In 



87r 



■In 



2tt 

TT 

+ 2 



A 



(41) 



p 



The form of the first term Fq in free energy H41II is similar to H29|l for an incompressible surface with 7 = 0. The 
second term Fi is the leading linear correction in the small variable eA/A*. The harmonic (Schulman) elastic energy 
contributes only to the quadratic correction. 

Our result (|41|l should be compared with the free energy H32|) of an incompressible surface in the limit 7 — s- 0. As 
has been discussed in the literature 3] (see also our discussion in section and as evident from our Eq. (|32|l , the 
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vanishing of 7 should not be confused with the vanishing of surface tension a which occurs when Ap (as we 

show in section HV CII . In the derivation of Ea. (|32|l . which apphes for the case of an incompressible surface and for 
which the total number of surface molecules is not fixed, we have used relation (PO)), which enforces the conservation 
of mass of the surface patches. Here, we consider a different scenario, where the total number of molecules is fixed. 
Consequently, relation (|30() should be replaced with Ea. (|14|l in which we fix iV, the number of patches dividing the 
surface. In order to calculate a by differentiating F with respect to Ap, it is better to consider the explicit expression 
for Fo and to perform the integral in Eq.|0TJ. Setting the lower limit Aq = we get 



47r 



In 



2^ 



-NkgT 



In 



(42) 



Apart from the constant —2NkBT which can be discarded (if temperature is fixed), the free energy depicted by the 
above equation is reminiscent of the free energy F = — JVfcsTln (^Ap/X^N^ of a two-dimensional ideal gas consisting 
of N particles confined in an area Ap. The only difference between the free energies of those two systems is that the 
effective de Broglie thermal wavelength in Ea. H42l) is A* = A-^Stt/Jk. While A is a constant that depends only on 
the mean number of molecules included in each patch (which is not affected by changes in Ap if we use Eq. (|14|l with 
a fixed value of iV), the effective thermal wavelength A* may be area-dependent. The origin of this dependence is 
the bending modulus k appearing in the expression for A*. In the quadratic approximation of the Canham-Helfrich 
Hamiltonian [Eq.JSJ], the elastic bending energy is expanded around the flat reference state of area Ap and, in general, 
K depends on the properties of the reference state. As the number of surface molecules is fixed, the variation of Ap 
leads to the variation of the area density of the reference state and, presumably, to the variation of n. The exact 
calculation of k requires the knowledge of the molecular interactions between the molecules [23, |36| , and is beyond 
the scope of this paper. We will, therefore, restrict our discussion in the reminder of this paper to surfaces for which 
the variation of k (over the range of Ap of interest) is negligible. In such a case. 



The first correction to this result is 



CTl 



dFi 
dA„ 



dFi 
deA 



In 



N 



+ 1 - 87r/3K 



Svr V 71- 



keT 
Al 



(43) 



(44) 



In the thermodynamic limit this is merely a logarithmic size correction to cto , and we thus conclude that for Ap ^ A* 



NkgT 



(45) 



The negative surface tension implies that the surface opposes its contraction. This result is expected since reduction 
of Ap requires the surface to bend more and as a consequence pay a higher energy cost in order to attain the Schulman 
area. Negative surface tension with the same origin was previously explained in the context of vesicles in Ref. (I5.j , and 
has been also observed in simulations of bilayer membranes [23, Issl l39l | . While the very fact that the surface tension 
becomes negative is not surprising, its magnitude (fcsT per microscopic unit area, i.e., of the order of I-IO ergs/cm^), 
as predicted by Ea. (|45|l . is strikingly large. Such a large and negative surface tension implies that undulations with 
wavelengths longer than 2TTy/ f3K,{Ap/N) should be unstable, because the coefficients of \hq\^ in the free energy become 
negative. This undulation instability will not show up in computer simulations where the size of the sample is small 
[37, 38, 39]. For real physical systems this highly compressed regime is unattainable. It can be preempted by a 
fluid-solid phase transition ^3| (in which case the shear modulus of the solid surface may act against the increase in 
the total area), or by a transition of molecules from the surface to the solution It should be noted here that a 
large negative surface tension of the same order of magnitude is also predicted by the incompressible surface model 
[see Ea. (|32|l . with 7 = 0]. However, it is probably worthwhile reemphasizing that surfaces are nearly-tensionless for 
cr ~ rather than for 7^0. 



B. Tense surfaces 



Consider the case when Ap > As- In this regime the saturated area is not attainable because A> Ap. For 
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(a condition which for ^/^Kacu, ^ 1 excludes only a small regime of Ap close to As), i.e. for AAs ^ 
—As {2f3KAas)^^^^ , the optimal area is very close to Ap so that Ea. (|35|l can be approximated by 

The above approximation is obtained by neglecting the exponent in the numerator of the expression in the central 
part of Ea. (|35|l . and expanding the exponent in the denominator close to A A ~ 0. We also omit the term AA on the 
r.h.s of Eq.(|2SI). The solution of Eq.(|^ is 

AA* = -^^^ ^ ^ ^ . (48) 

SttPKaAAs 2PKaAAs 2f3KAasAAs ^ ' 

(Note that for the range of A As considered here the solution satisfies A A* ^ |AAs|, and so our approximation of 
the r.h.s of Eq.jSHl) is justified.). Combining Eqs.Hnj, (EEl, (ESJ, JIZI), and (gHl), and recalling that AA* < \AAs\, we 
find that the free energy is given by 



F ^ ^(AAs) 



knT- — ^ [ — In 



{(3{KA\AAs\/As)q' + f3Kq'') PX^ 
2-K 



(49) 



If temperature is fixed then the second term in the above equation is a constant and can be ignored. The first term 
depicts the energetic contribution of Schulman elasticity to F . This is simply the elastic energy of a flat surface having 
been stretched to an area Ap > As- The surface tension experienced by such a fiat surface is 

ao = -^AAs = ^\AAs\ > 0. (50) 

The third term in Ea. H49(l represents the entropic part of F. It is similar to the second term in Ea. H29|l . which is the 
leading thermal correction to the free energy of an incompressible surface, with 7 = erg. The (full) derivative of this 
term with respect to Ap yields the correction to the linear relation ct = ctq, between the surface tension and A As- It 
can be shown that for Kao-s/ 1^ 2> 1 [see relation (|34() ]. the sign of this correction is positive which means that thermal 
fluctuations lead to a super-Hookean elasticity. This observation can be understood by noting [see Ea. H48|l ] that the 
optimal area A* is independent of k, and that it gets closer to Ap by increasing |AAs|- The approach of A* to Ap 
involves a strong suppression of the thermal fluctuations, and that naturally reduces the entropy of the surface and 
thereby increases the surface tension. As already said, this entropic surface tension is also found for incompressible 

surfaces [n 113, mill mm. Stretching experiments of vesicles are often analyzed assuming the decoupling of 

the elastic and entropic surface tensions. However, our Ea. l|l^ suggests the existence of a rather complicated interplay 
between the elastic and the entropic contributions, which is expressed by the fact that the stretching modulus Ka 
appears in the integral expression for the entropic component. 



C. Schulman surfaces 



Finally, we discuss the case when Ap ~ As. More precisely, we consider the regime where the projected and 
saturated areas are sufficiently close to each other that 



V 2(jKAas 



where A As is either positive or negative. Instead of Ea. (|47|l we now have 



/\2 /I 



Up to second order in the parameter AAs/lvDAs) < 1, the solution of this equation is 
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and hence 



Us (AA*) = pKj 



D 



^ 2As ?>y/DAl^ 

Eas. p9|l . (|26|) . H53|) . and H54|) yield the foUowing form for the free energy 



F ~ ^{AA^f ^ DKaAs 



(54) 



(55) 



Ap rdq 

(27r)2Ao 2 



In 



/3Z'A2g2 ( ^^^ KaAA, , if^(AA 



27r 



2A, 



The last three terms in the argument of the logarithmic function are smaller than the first one [see relations H34|) and 
(|51|l ]. Using a Taylor expansion of the logarithm and keeping only terms up to first order in Pk/^/PKaAs and second 
order in AAs/{\/DAs) we obtain 



F o^^{AA,f-DKAAs + kBT ^ 



dq 

(2^)2 ./n T 



(2^)^7o 2 
-KaAA, , ifA(A^ 



2A,, 



27r 
{KaAAs 

^1^ 



2 



KaVD 



-l + \n(2pKAX^VD 



k,T 



(2^)2 



dq 



-KaAA, , ifA(A^ 



2A,, 



9,^/DAl 



KaVD 



(KaAAs 



K,q 



KaVD 



(56) 



If we fix the temperature and ignore the small dependence of k on Ap, then the second and third terms in Ea. l|5t)|) 
are constants which can be dismissed. The first term is the zero-temperature Schulman elastic energy, while the forth 
term is the thermal correction. The integrals in the latter can be easily performed, and after several mathematical 
manipulations we find that (up to an irrelevant additive constant) 



Ka 
2As 



(AA^f 1 



IGttkVD 

KaUs 



KaVdAAs 



KAas 



(57) 



In the A As regime which we discuss in this section [see Eq.(|5TJ] the second term in Ea. (|57() should not be considered 
as correction to the first, but rather as the dominant term in the free energy. From Ea. H57|) we conclude that thermal 
fluctuations leads to a shift in the value of the projected area at which the free energy is minimal. The free energy 
attains its minimum at 



a,, = AsIi-Vd 



o 



Kaus 



A 



off 



(58) 



which can be regarded as the effective Schulman area at which a ~ 0. The very fact that Af^ < Ag is not surprising 
since the total area is always larger than the projected area. For typical values of phospholipids: k = lOfc^T ^ 5x10^^'^ 
ergs, Ka ~ 200 ergs/cm^, and ~ 10~^^ cm^ (the area of a patch consisting of 1-2 molecules), we find that ^/D ~ 0.1 
and so the shift in the equilibrium projected area is about 5% of Ag. We also conclude from Ea. H57|) that thermal 
fluctuations increase the compression-stretching modulus Ka- The effective modulus is 



K 



cff 



Ka 



IGttkv D 



(59) 



For the values of the relevant quantities quoted above (characteristic of phospholipid bilayers) we find, quite surpris- 
ingly, that the magnitude of the thermal stretching modulus is comparable to the bare modulus Ka- 
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In the Evans- Rawicz experiment ^J, linear elastic response was measured for areal strain rj = (Ap — ^4^^^)/^°*^ < 
0.05. This is precisely the regime studied in the present section, and so our theoretical discussion explains well the 
experimental behavior. For extremely small values of the (< 0.005) the surface tension grows exponentially with 
the strain. This elastic behavior has been recently explained in the framework of the incompressible surface model, 
and has been attributed to the dependence of the entropic surface tension on the optically visible area (2J| . 

V. SUMMARY 

We have studied the elasticity of fluctuating two-dimensional compressible (elastic) surfaces. Starting with a 
Hamiltonian including the Canham-Helfrich bending energy and an area-dependent term, we have first developed a 
general formalism for analyzing the statistical-mechanical behavior of such systems. We have used the formalism to 
reproduce the expression for the temperature-dependent tension of incompressible surfaces. Our investigation of the 
elastic behavior of compressible surfaces reveals that the minimum of the free energy is obtained when the projected 
(frame) area Ap = Af^ , where A'f^ is slightly smaller than the saturated area Ag. In the vicinity of A'f^ the elastic 
energy is depicted by the Schulman quadratic form with an area elasticity modulus K'^ slightly larger than the 
bare modulus Ka- The response of the surface to compression and stretching away from Af^ is very different. The 
stretching behavior is dominated by the area elasticity. Thermal fluctuations introduce a correction to Hooke's law, 
which increases the surface tension. By contrast, area elasticity plays a very small role in the compression behavior 
because the total area of the fluctuating surface is (effectively) constrained to Ag. The negative surface tension 
observed in this regime should be attributed to fact that upon compression, the total area stored by the fluctuations 
(the excess area) grows. The enhancement of thermal fluctuations leads to an increase in the bending energy of the 
surface. 

We conclude by briefly reviewing the factors which have been left out of our model, and which influence the elastic 
behavior of the surface. One of them is the exchange of molecules between the surface and the solution. It can be 
dealt with within the framework of a grand-canonical ensemble where both the area A and the number of patches N 
are allowed to fluctuate. Understanding the dependence of the bending modulus k on the density of the reference state 
(i.e., on the ratio Ap/Ag) is another challenge. This is a complicated matter since the values of the phenomenological 
parameters k, Ag, and Ka are determined by the structure of the molecules and the interactions between them. It 
would therefore require a theoretical study at the molecular level (see, e.g. ^3 HE 113 )■ For strongly fluctuating 
surfaces we should include corrections to the quadratic form of the Canham-Helfrich Hamiltonian and to go beyond the 
lowest order in ksT/n. In that case "measure factors" 21j which rather complicate the mathematical treatment must 
be introduced in order to correct the partition function. Finally, we mention the possibility that strongly-stretched 
membranes can reduce their elastic energy by developing pores. A discussion at this subject can be found elsewhere 

m El mill El m 
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